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Abstract 

A mathematically and physically sound three- 
degree-of-freedom dynamical model that emu- 
lates low- to high-confinement mode (L-H) tran- 
sitions is elicited from a singularity theory cri- 
tique of earlier fragile models. We construct a 
smooth map of the parameter space that is con- 
sistent both with the requirements of singular- 
ity theory and with the physics of the process. 
The model is found to contain two codimension 
2 organizing centers and two Hopf bifurcations, 
which underlie dynamical behavior that has been 
observed around L-H transitions but not mir- 
rored in previous models. The smooth traver- 
sal of parameter space provided by this analysis 
gives qualitative guidelines for controlling access 
to H-mode and oscillatory regimes. 

I. INTRODUCTION 

A unified, low-dimensional description of the 
dynamics of L-H transitions Q would be a valu- 
able aid for the predictive design and control of 
confinement states in fusion plasmas. In this 
work we report significant progress made toward 
this goal by developing the singularity theory ap- 
proach to modeling L-H transitions that was in- 
troduced in Q. The results give new insights 
into the role of energy exchange and dissipa- 
tion in the onset, evanescence, and extinction of 
discontinuous and oscillatory action in confined 
plasmas. 



The title of this paper refers to the philoso- 
phy of singularity theory || as applied to dy- 
namical models: that paths through parame- 
ter space should be smooth and continuous, and 
that parameters should be independent and not 
fewer than the codimension* of the system. 

Since 1988 Q many efforts have been made 
to derive unified low-dimensional dynamical 
models that mimic L-H transitions and/or as- 
sociated oscillatory behavior [§-21]. All of these 
models have contributed to the current view in 
which the coupled evolution of poloidal shear 
flow and turbulence establishes a transport bar- 
rier. However, as was shown in [Q, the mod- 
els often founder at singularities. Consequently, 
much of the discussion in the literature concern- 
ing the bifurcation properties of L-H transition 
models is qualitatively wrong. 

We examine the bifurcation structure of a 
semi-empirical dynamical model for L-H tran- 
sitions ||, and find it needs two major opera- 
tions to give it mathematical consistency: (1) 
a degenerate singularity is identified and un- 
folded, (2) the dynamical state space is expanded 
to three dimensions. 

We then analyse the bifurcation structure of 
the enhanced model obtained from these opera- 
tions, the BD model, and find it consistent with 
many known features of L-H transitions. In par- 
ticular, this is the first model that can emulate 
the onset and abatement of oscillations in H- 
mode, and direct jumps to oscillatory H-mode 



|22,|3|. 
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II. BIFURCATION STRUCTURE OF 
THE DLCT MODEL 

This paradigmatic 2-dimensional model || 
comprises an evolution equation for the turbu- 
lence coupled with an equation for the flow shear 
dynamics derived from the poloidal momentum 
balance: 



that F oc v' 2 . 
is 



Thus, the modified DLCT model 



— = 7 N - aFN - (3N 2 
at 

dF 

— = a FN - uF. 
dt r 



N is the normalized level of density fluctuations, 
F is the square of the averaged ExB poloidal 
flow shear. The fluctuations grow linearly with 
coefficient 7 and are damped quadratically with 
coefficient [3. The exchange coefficient a is re- 
lated to the Reynolds stress, and the damping 
rate \xF is due to viscosity. ^ 

Following the procedure outlined in [Q] we 
form the bifurcation function g = X(F,~f), and 
identify the singular points where g = gF = 
0. We find the unique physical singularity 
(F, r y)T = (0, (3(j,/a), which satisfies the addi- 
tional defining conditions for a transcritical bi- 
furcation: 



9j = 0, Qff + 0, det d 2 g < 0, 



(3) 



where det d 2 g is the Hessian matrix of second 
partial derivatives with respect to F and 7. 
Evaluating (|3|) at T gives g 7 = 0, gpF = 
—2a 2 /(3, det d 2 g = -a 2 / (3 2 . The bifurcation di- 
agram showing the transcritical point T is plot- 
ted in Fig. |l]a. (In this and subsequent diagrams 
stable solutions are indicated by continuous lines 
and unstable solutions by dashed lines.) 

However, Fig. [T]a does not represent the com- 
plete bifurcation structure of the DLCT model 
because of the following generic property of the 
transcritical bifurcation: it is non-persistent to 
an arbitrarily small perturbation. Since the 
poloidal shear flow v' is symmetric under the 
transformation v' — > —v' it is appropriate to 
introduce the perturbation term pF 1 ! 2 . (Note 



dN , 
— = 7 AT - aFN - (3N 2 



dF 
~dt 



aFN - {iF + ipF 



1/2 



(1) 
(4) 



The perturbation term in Eq. |3] represents 
a physically inevitable source of shear flow that 
breaks the symmetry of the internal shear flow 
generation and loss rates. The physics comes 
from non-ambipolar ion orbit losses that produce 
a driving torque for the shear flow. This con- 
tribution to the total shear flow evolution can 
be quite large p^fl , and in fact the early models 
for L-H transitions relied exclusively on a large, 
nonlinear ion orbit loss rate [|],||]. However, ion 
orbit loss alone cannot explain turbulence sup- 
pression. Here we treat this term as part of a 
more complete picture of L-H transition dynam- 
ics, and emphasize its symmetry-breaking na- 
ture by assigning the simplest consistent form to 
it while recognizing that ip may be a nonlinear 
function 97(C), where £ may include dynamical 
variables and parameters. 

Some bifurcation diagrams for increasing val- 
ues of <p are plotted in Fig. [l]b-d. We see im- 
mediately that solution of one problem causes 
another: a nonzero perturbation term does in- 
deed unfold the degenerate singularity T, but it 
releases another degenerate singularity T s . 

Before proceeding with a treatment of the 
new bifurcation T s we highlight three important 
issues: 

1. Since 99 is inevitably nonzero in experiments, 
no transition can occur at all in the vicinity of 
T, neither first-order or second-order, contrary 
to what is stated in ||. 

2. Both N and F change continuously with 7 in 
the same direction. The fact that T is a trans- 
critical bifurcation tells us that only two param- 
eters — (p and any one of the other parameters 

- are required to define the qualitative struc- 
ture of the problem. The bifurcation diagram 
with N as state variable is plotted in Fig. 0, 
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which should be compared with Fig. ||c. As it 
stands, the model therefore cannot emulate tur- 
bulence stabilization by the shear flow, contrary 
to what is stated in ||. 

3. To ascertain whether the model can exhibit 
periodic dynamics as stated in || we look for 
a pair of purely complex conjugate eigenvalues. 
For Eqs |l] and ^ the defining conditions for Hopf 
bifurcations may be expressed as 



g = tr J = 0, det J > 0, 



d_ 



tr J ^ 0, (5) 



where J is the Jacobian matrix. We find that 
det J < where the equalities in Eq. || are ful- 
filled, therefore oscillatory dynamics arising from 
Hopf bifurcations cannot occur. 

This does not rule out the possible existence 
of periodic behavior arising from rare and patho- 
logical causes. According to Dulac's criterion 
[25] Eqs [j] and [I] possess no periodic solutions 
arising from any cause if there exists D such that 
the quantity 



(6) 



never changes sign. Here W = W(N, F) = 
dN/dt, Y = Y(N, F) = dF/dt, and the Dulac 
function D = D(N,F) is a real positive func- 
tion. Choosing 3 = 1 we find that 



S = a(N — F) — 2N/3 + 7 - \i + ^F' 1 / 2 /2, 



(7) 



which clearly can switch sign. However, there 
may exist a more exotic Dulac function that for- 
bids a change of the sign of S. We have not found 
oscillatory solutions numerically in this system. 

Returning to the new singularity T s we 
find that it is also a transcritical bifurcation: 
the conditions (||) evaluated at (F, 7)^ = 
(<p 2 /fi 2 , a>ip 2 /fj, 2 ) yield g 1 = 0, g FF = 

-/3^ 6 /(2aV 4 )-^ 3 /^ 2 , det A = -^7(4«V)- 
Does the DLCT model therefore require a second 
perturbation term, this time to Eq. |l], to unfold 



We remark here that often there is more than 
one universal unfolding for a given bifurcation 
problem, and we turn to the physics to decide 
which is physically consistent. For the pertur- 
bation in Eq. || that unfolded T we chose the 
form ipF 1 / 2 because it is physically inevitable 
that the symmetry v' — > —v' be broken. How- 
ever, there is no matching physics for a simi- 
lar term in Eq. |l]. Another possibility is that 
T s is spurious, created by an unwarranted col- 
lapse of a larger state space. This idea leads to 
a suggestion that is supported by the physics, 
that another dynamical variable is intrinsic to 
a low-dimensional description of L-H transition 
dynamics. 

III. INTRINSIC 3-DIMENSIONAL 
DYNAMICS OF L H TRANSITIONS 

We introduce the the third dynamical vari- 
able by assuming that 7 = j(P), where P is 
the pressure gradient, as have a number of other 
authors ]p|, p^ , ^3| , p^j20[| . Assuming the simplest 
evolution of P and that j(P) = jP, we arrive at 
the following augmented model, obtained purely 
from dynamical and physical considerations: 

>lf q-^PN 



dt 
dt 



aFN - f3N 2 
- LiF + ipF 1 / 2 . 



(8) 
0) 



In Eq. |H q is the power input and e is a di- 
mensionless parameter that regulates the contri- 
bution of the pressure gradient dynamics to the 
overall evolution. The dynamics is essentially 3- 
dimensional with e ~ O(l), but for e <C 1 or 
e S> 1 the system can evolve in two timescales: 
1. The original "slow" time t. For s — > 0, 
edP/dt ~ and P ~ q/{^N). The system col- 
lapses smoothly to 
dN 

~dT 
dF_ 

~dt 



q - aFN - /3N 2 
aFN - fiF + LpF 1/2 . 



(10) 
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The organizing center is the unique transcritical 
bifurcation (F,q,tp)T = (0, /3/x 2 /a 2 , 0), the spu- 
rious T s is non-existent, and there are no Hopf 
bifurcations. For e > 1 we define 6 = 1/e and 
multiply Eq. |§] through by 5; taking the limit as 
5 — > gives dP/dt m 0, from which P = Pq. We 
recover the same form as Eqs |l| and |], 



dN 

~dt 
dF 

~d£ 



>yP N - aFN - f3N 2 
aFN — fiF + (fF 1 / 2 , 



(!') 



along with the "good" bifurcation T and the 
"bad" bifurcation T s - - therefore we suggest 
that this is a non-physical limit for e. 
2. In "fast" time r = ejt and, recasting the 
system accordingly, it can be seen that on this 
timescale the dynamics becomes 1-dimensional 
in P in both limits. 

The organizing center of the bifurcation 
problem obtained from Eqs H ||, and || is 
the unique transcritical bifurcation (F, q,(p)r = 
(0, Pn 2 /a 2 , 0), g FF = -a 2 / (3, det d 2 g = 
— a /(4/? 2 // 2 ), and the spurious singularity T s is 
non-existent. We now have the bones of an im- 
proved dynamical model for L-H transitions, but 
it still does not emulate the following character- 
istics of L-H transitions: (a) Hysteresis: Since 
there is no non-trivial point where gpF = it 
cannot model discontinuous transitions or clas- 
sical hysteretic behavior, (b) Oscillations in H- 
mode: These have not been found numerically. 
In a 3 - dimensional dynamical system it is, of 
course, very difficult to prove that oscillatory so- 
lutions do not exist. 

Evidently we need more nonlinearity or 
higher order nonlinearity to produce enough 
competitive interaction. To obtain multiple 
solutions, at least, the bifurcation equation g 
should map to the normal form for the pitchfork 
bifurcation h = ±x 3 ± Ax. 

Several authors have taken the viscosity co- 
efficient ELS cl function of the pressure gradient, 
but usually it is treated constant. In 
the viscosity was considered to be the sum of 



neoclassical and anomalous or turbulent contri- 
butions, both with separate power-law depen- 
dences on the pressure gradient. We shall adopt 
this bipartite form and in Eq. take 



fJL = MP) = fineoP n + ^nP" 



(11) 



Equations §, | and | with flO]) comprise the BD 
model. 

The values of the exponents n and m are not 
precisely known empirically or from theory. In 
[26] fj, an is given as having a P 3 / 2 dependence, 
but is also subject to the additional influence of 
a P-dependent curvature factor. In this work we 
take n = —3/2 as in Jl^] and m = 5/2. 



IV. BIFURCATION STRUCTURE OF 
THE BD MODEL 

The bifurcation problem obtained from the 
BD model contains two codimension 2 organiz- 
ing centers: 

1. The defining conditions for the pitchfork, 

9 = 9F = 9FF = 9q = 0, 9FFF + 0, g F q / 0, 

(12) 

find this singularity occurring at (F,q,/3,cp) p = 
(0, 8^ 8 /in / co7/(7 7/8 a), (7 3 / 8 a 7 )/(8^Vn / eo), 0), 

9FFF = -12(7 7 / 8 )/iI^ 8 / Uneo7/a> 9Fq = 

2(7/Li an /^ ne o) 1//4 . The pitchfork p becomes a 
transcritical bifurcation T l away from the criti- 
cal value of (3. 

2. Another transcritical bifurcation T u occurs 
at (F,q,ip) T u (0, P 2 7 //?, 0), 
9ff = -2P 7 3 (7P 5 / 2 a 7 - 8^ nco /3)/(P 7 / 2 Q/3 7 2 ), 
det d 2 g = -(-3P 5 / 2 q 7 + 8/i neo ^) 2 /(4P 7 Q 2 7 2 ). 

T l and T u are annihilated at a second codi- 
mension 2 bifurcation. The defining conditions 
for this point are 

9 = 9F = 9q = det d 2 g = 0, g F F + 0, g Fq + 

(13) 

At this point we find 



4 



(0,(8(7 1 / 8 )/i nC o7)/(3a(^nco/Aian) 1/8 ), 

3(^neo/Man) 5/8 a7/(8(7 5 /Vnco)),0), 

9ff = -64(7/i an ) 5 / 8 /i^o7/(3a), gF g = 
4(7^ an /^ neo ) 1 /4. 

In Fig ||a the partially perturbed bifurcation 
diagram is plotted, showing the lower and up- 
per transcritical bifurcations T l and T u . In Fig. 
||b the fully perturbed, physical bifurcation dia- 
gram is plotted, where 92 > 0. There are also two 
Hopf bifurcations on the upper H-mode branch 
in Fig. H linked by a branch of stable limit cycles. 
The dotted lines mark the maximum and mini- 
mum amplitude trace of the limit-cycle branch. 
This reflects the passage through an oscillatory 
regime that is often observed in experiments. 

Since it is a codimension 2 bifurcation prob- 
lem, the qualitative structure is fully defined by 
q and two auxiliary parameters. One of these 
is obviously ip, the other may be any one of the 
other parameters. We choose (3 because we are 
interested in the effects of poor turbulence dis- 
sipation (i.e. low /?). Figure |] illustrates how a 
jump can occur directly to oscillatory states, a 
phenomenon which is frequently observed. 

Figure [|, to be compared with Fig. ||]b, shows 
that the BD model does indeed reflect shear flow 
suppression of turbulence. 

V. DISCUSSION AND CONCLUSIONS 

A dynamical model that emulates much of 
the typical behavior around L-H transitions has 
been elicited from an earlier fragile model that 
had serious flaws by considering the relationship 
between bifurcation structure and the physics of 
the process. Built in to this model are the follow- 
ing major dynamical features of L-H transitions: 

1. Discontinuous, hysteretic transitions, or 
smooth changes with power input, depending on 
the degree of turbulence dissipation 0, or equiv- 
alently, the viscosity. 

2. Two Hopf bifurcations in H-mode. It is the 
first model that can emulate the onset and abate- 
ment of oscillatory behavior, and a transition di- 



rectly into oscillatory H-mode. 

3. Turbulence suppression by the shear flow. 

4. A maximum in the shear flow generated by 
the turbulence, followed by a decrease as the 
power input flowing to the turbulence is raised. 

5. Turbulence generation from non-ambipolar 
losses. 

Finally, we note that the existence of two 
codimension 2 bifurcations is suggestive: Should 
there be an expansion of the system, perhaps ex- 
pressing fluctuations of the magnetic field, that 
creates (or annihilates) the two bifurcations at a 
codimension 3 singularity? In other words, does 
a more complete model contain an organizing 
center of higher order? In singularity theory we 
persevere in seeking higher order behavior: that 
is how the relationship between a model and the 
process it represents is tracked. This question is 
currently under investigation. 
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FIGURES 





FIG. 1. Bifurcation diagrams of the DLCT model showing how a perturbation of the shear 
flow unfolds T but introduces T s . a = 1, (3 = 0.77, \i = 1. a. y> = 0, b. (p = 0.05, c. <p = 0.5, 
d. 99 = 1. The region F < is not within the phase space but is included to make the nature of T 
clearer. 



7 



0.8 r 




✓ 



A I I I I 

0.5 1 1.5 2 

n Y 

FIG. 2. Same as Fig. [y(c) except with N as the state variable. The region N < is not 
within the phase space but is included to make the nature of T s clearer. 
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FIG. 3. Bifurcation diagrams BD model, a = 2.4, (3 = 1, 7 = 1, e = 1, /i neo = 1, 
A*an = 0.05,n = —1.5, m = 2.5. a. Bifurcation structure of the partially perturbed system, 
with if = 0. b. <£> = 0.05. For clarity the lower unstable branches are not plotted in (b). 
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FIG. 4. Under extreme conditions, where the turbulence dissipation rate is low relative to the 
generation rate, the jump at the lower limit point can occur directly to an oscillatory state on the 
H-mode branch. (3 = 0.1, other parameters as for Fig. 0. 
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